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1 .  Introduction 

The  ways  In  which  molecular  vibrations  relax  at  surfaces  are  of  considerable 

interest  in  dynamical  studies,  not  least  for  investigating  the  possibility  of 

bond-selective  chemistry.  Recently,1  an  experiment  was  reported  which  has  for 

the  first  time  directly  measured  the  rate  of  vibrational  energy  relaxation  of  an 

adsorbed  species  (OH  on  Si02,  v-1+v-O).  Knowledge  of  the  total  width  of  the 
2 

absorption  line  then  allows  an  estimate  of  the  pure  dephasing  rate  through  the 
relation  1 /T2  -  2/T1  +  1/T^,  where  T^  and  are  the  energy  and  pure  phase 
relaxation  times,  respectively.  There  are  undoubtedly  many  relaxation  pathways, 
including  direct  coupling  of  the  molecular  mode  to  the  phonons  and  formation  of 

3 

electron-hole  pairs.  Silica,  however,  is  a  substrate  in  which  the  latter 

mechanism  can  be  effectively  ruled  out. 

We  advance  here  a  treatment  of  pure  vibrational  dephasing  of  an  adsorbed 

diatomic  molecule,  via  direct  vibration-phonon  coupling.  For  the  sake  of 

comparison,  we  give  an  expression  for  a  T1  pathway  appropriate  for  0H/Si02.  The 

model  involves  an  anharmonic  oscillator  (A-B)  coupled  through  an  anharmonic  bond 

(S-A)  to  a  linear  chain  of  harmonic  oscillators.  We  show  that  phase  relaxation 
2 

is  about  10-10  times  faster  for  this  sy stern  than  energy  relaxation,  in  broad 
agreement  with  experiment.  We  also  find  the  coupling  to  low  frequency 
librational  modes  to  be  unnecessary  for  explaining  experimental  results.  In 
this  regard  we  feel  it  worth  emphasizing  the  advantages  of  a  linear  chain  model 
when  the  admolecule  is  aligned  end-on  at  the  surface.  The  assumption  is  that  the 
dynamics  of  admolecule  and  surface  is  dominated  by  the  motion  of  the  whole 
system  in  the  z-direction.  Such  a  model  has  recently  been  successfully  applied 
to  the  study  of  the  desorption  of  an  adatom.  The  incorrect  notion,  that  a 
linear  chain  model  is  one-dimensional  in  every  respect,  is  quite  widespread.  In 
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fact,  the  density  of  states,  a  quantity  needed  for  the  calculation  of  rate 
information,  is  obtained  from  a  fully  three-dimensional  (Debye)  treatment,  in 
which  the  phonon  modes  are  assumed  to  form  an  isotropic  continuum.  The  only 
real  assumption  apart  from  this  is  that  the  molecular  dynamics  is  dominated  by 
motion  in  a  particular  direction.  Provided  this  is  the  case,  such  models  are 
expected  to  give  the  correct  order  of  magnitude. 

2.  Model 

The  details  of  the  model  are  as  follows.  The  momentum  and  position  normal 

coordinates  of  the  diatom-surface  system  are  denoted  by  (p^q^)  and  (p2,q2). 

When  the  frequencies  of  the  bonds  are  sufficiently  different,  the  motion  of  the 

center  of  mass  of  AB  relative  to  the  surface  may  be  identified  as  one  normal 

mode,  and  the  vibration  of  A  relative  to  B  as  the  other.  Thu3  we  may  make  the 

identities  q1  -  z  and  q2  -  x.  We  shall  focu3  on  these  position  coordinates  for 

the  Taylor  expansion  of  the  potential,  although  in  fact  the  nonadiabatlc  terms 

5 

due  to  p1  and  p2  give  contributions  of  a  similar  order  of  magnitude.  The  first 

term,  in  the  Taylor  expansion  of  the  interaction  potential,  which  gives  rise  to 

pure  dephasing  (i.e.,  the  leading  term  in  the  expansion  which  is  diagonal  in  the 

2  6 

vibrational  mode  of  the  admolecule)  is  the  quartlc  k22x  term. 

In  the  theory  of  spectral  lineshape,  the  width  of  a  line  is  given  by  the 

7 

imaginary  part  of  the  Zwanzlg  memory  tetradic  <T(uf,i  +  ie)>fi  where  c  is  an 

infinitesimal  quantity  and  the  transition  is  being  made  between  the  initial  (i) 
and  final  (f)  state  of  the  system.  This  may  be  expanded  in  powers  of  the 
interaction  Liouville  operator.  In  the  Markov  limit  (short  memory,  T  is 
constant)  and  in  first  order,  there  are  two  kinds  of  term:  one  represents 
contributions  to  the  linewidth  from  inelastic  scattering  out  of  levels  f  and  i, 
the  other  represents  elastic  transitions.  The  latter  are  identified  with  the 

g 

"pure  dephasing"  interaction  and  an  associated  rate  1 /T2  given  by: 


k-  1  !,  <‘o6<l<«l*“>Vla  -  <v(b)>2V  «(-„,)  . 

2  “fi  a,  B 


where  the  bath-system  interaction  is  assumed  to  be  of  the  Form  in 

obvious  notation. 


fi  m 


is  the  relative  displacement  of  the  levels  i  and  f  of  the  molecule  due  to  the 
anharmonic  interaction  with  the  bath;  d  is  the  canonical  density  operator  for 
the  bath,  which  is  assumed  to  be  diagonal  and  therefore  describes  the  bath  at 
thermal  equilibrium.  The  indices  u  and  B  represent  entire  states  of  the  bath. 
Put  another  way,  Eq.  (la)  represents  the  thermally  averaged  rate  of  Rayleigh 
scattering  of  phonons  by  the  system  (admolecule) . 

On  the  other  hand,  the  inelastic  rates  are  given,  in  tfte  Markov  limit,  by 


the  Golden  Rule  expression: 


l  |<a|V(bV|2  l4f  I2  • 


t, 


which  is  essentially  the  thermally  averaged  rate  of  phonon  absorption  and 
emission  by  the  system. 

In  order  to  Introduce  the  bath-system  coupling,  we  first  evaluate  the 
dynamics  of  AB  moving  with  respect  to  a  frozen  surface  atom  S  (i.e. ,  the  surface 
is  represented  by  an  atom  of  infinite  mass).  The  influence  of  the  thermal 
fluctuations  of  the  lattice  is  then  introduced  by  allowing  the  surface  to 
displace  by  an  amount  u,  which  is  fluctuating  with  time.  Moreover,  the 
assumption  of  Markovian  behavior  implies  u(t)  is  uncorrelated  with  u(t'),  t*t'. 
We  replace  z  with  z-u  and  express  u  in  terms  of  a  Fourier  sum  over  the  normal 


coordinates  of  the  lattice:" 


...  ,  •fi  .1/2.  -1/2  ,  t  "‘“k1 

ult)  -  I  « k  <V 

k 


iw  t 

+  V  )  * 


(3) 
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where  there  are  N  atoms  in  the  chain,  each  of  mass  M,  and  ak,ak  are  the  creation 

and  annihilation  operators  for  mode  k,  which  has  frequency  w  . 

First,  we  consider  the  dephasing  of  a  two-level  molecule  attached  to  a 

surface.  By  making  a  Taylor  expansion  of  the  interaction  it  is  possible  to  pick 

out  terms  which  contribute  to  the  process  of  interest.  The  fir3t  term  in  the 

2  2 

interaction  which  gives  rise  to  pure  dephasing  Is  k^  u  x  ,  and  therefore 

V(b)v(3)  =  k2Q2  uV  .  (4) 

Here,  the  subscripts  represent  the  degree  of  excitation  of  phonons,  adbond  and 
molecule,  respectively.  Substituting  Eqs.  (3)  and  (4)  into  Eq.  (la)  yields 
lk2 

Tj  "  ^2^2  k^k,d8i*  (l<aWak'  +  akak  •  I  6>  I  “  nkl5kk')  A1C>)<S((V<V)  *  (5a) 

—  9^1 

Here,  is  the  thermal  average  of  nk»  and  is  equal  to  (e  -1)  ,  with  0  - 

K'V' 

Converting  the  sum  to  an  integral  gives 


1  ™22  (s)f 

T  '  ‘  ~T1  A10  I  dukdu,k' 

2  4M  N  J 


»(uk)p(uk*)  - 

— — ; - nk(nk«+1)  • 

ukV  k  k 


where  p(x)  is  the  density  of  modes  of  frequency  x.  We  now  invoke  the  Debye 

approximation,  in  which  the  phonon  states  of„a  solid  form  an  isotropic  continuum 

3wk 

with  a  density  of  states  given  by  p(u»k)  -  — ~  N,  where  wD  is  the  Debye  frequency 

UD 

of  the  solid.  This  yields: 

!_  9vk22  (s)  (6  . 

T’  2  6  A10  1 2  ’  (6a) 

2  5M  Ujj 

where 

I  -  j“D  du^n*  ♦  nk)u,^  .  (6b) 

( s ) 

In  order  to  evaluate  in  Eq.  (6a),  we  take  the  states  |l>  and  |0>  to  be 

lu  1  m  1  m 

states  of  a  harmonic  oscillator,  when 


■*  "  <  *  * 
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<j|x2|j>  «  — -ft  (7) 

m  m  2viab“ab 

Here,  uAQ  i3  the  reduced  mass  of  the  admolecule  and  u>AB  Its  fundamental 
frequency.  Moreover,  we  can  define  two  limiting  cases  for  which  analytic 
solutions  to  Eq.  (6b)  are  possible.  In  the  limit  of  low  temperature,  0  >>  1 , 
and  u>D  may  be  taken  to  be  infinity.  In  the  limit  of  high  temperature,  e  <<  1 
and  the  exponentials  in  the  expression  for  n^  may  be  expanded,  while  retaining 
only  the  first  term.  Then  we  arrive  at  the  following  equations  for  the  pure 
dephasing  rate: 


o  3,2  ,3 

3lT  k02kB  3 

4fiM2p2  w2  J3 
A6  AG  D 


kBT  «  fuaD 


1_  9lf  k02kB  2 
T'  “  ,M2  2  2  5 

2  "n  ^AB^AB^D 


kBT  »  H 


Consider  now  the  energy  relaxation  of  A-B.  In  the  case  of  OH/SiO^,  u>D  -  3^0 

cm  1 ,  uAS  »  970  cm  1  and  uAB  -  3700  cm  1  .  Therefore,  the  first  important 

5 

interaction  giving  rise  to  energy  relaxation  is  the  sextic  term,  (z-u)  x,  and  in 
4 

particular  uz  x  which  destroys  one  phonon,  creates  four  quanta  of  adbond 
excitation  and  destroys  the  quantum  of  vibrational  excitation  in  OH.  Moreover 
we  note  that  a  direct  one  phonon  deactivation  of  AB  is  not  allowed  by  energy 
conservation,  and  that  it  would  require  at  least  eleven  phonons. 

We  now  identify  eigenstates  of  the  system  as  products  of  adbond  states  and 
molecular  states.  Using  Eq.  (1b)  together  with  the  Debye  density  of  states 
y ields 


r~  “  n(w  )w  0(o)  "WfJ  , 

ll  P  ”  P  0 


(9a) 
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where 

T  •  kU1  a<lll2Vam«>M,>.  .  (9b) 

1) 

and  uip  ■  ^“sa  -  UAB;  kl4l  is  the  coe^^icient  of  112  x>  For  simplicity  we  have 
taken  the  eigenstates  of  the  adbond  to  be  harmonic  oscillator  states.  This  is 
clearly  a  cruder  approximation  than  making  AB  harmonic,  especially  given  the 
level  of  excitation  of  SA  (four  quanta).  However,  there  is  little  detailed 
Information  about  the  silica-hydroxyl  bond,  and  we  are  in  any  case  interested 
only  in  making  an  estimate  of  the  rate  of  energy  relaxation  for  this  system. 

The  introduction  of  anharmonicities  into  the  SA  bond  will  make  the  calculation 
only  slightly  more  complicated.  (Clearly,  the  interaction  Hamiltonian  depends 
only  on  the  anharmonlcity  of  the  bond!) 

It  remains  to  determine  k and  k^.  Obviously,  we  must  here  take  a 
proper  account  of  the  anharmonlcity  of  the  adbond,  SA;  indeed  we  shall  make  the 
assumption  that  this  is  the  principal  contributor  to  the  coupling  between  the 
adbond  and  the  molecular  band  AB.  The  validity  of  this  assumption  will  be 
assessed  at  a  later  time.^  Further  we  make  the  assumption  that  z,  the  distance 
between  the  surface  and  the  center-of-mass  of  AB,is  given  approximately  by  the 
distance  SA.  This  assumption  is  expected  to’  be  reasonable  when  A  is  must  more 
massive  than  B  (as  in  OH).  Therefore,  the  motion  of  the  AB  bond  gives  rise  to  a 
modified  Morse  potential  for  SA  in  which  the  distance  parameter  is  z-ox,  where  o 
is  the  ratio  A  simple  expansion  of  the  Morse  potential  around  the 

equilibrium  distance  then  gives  terms  containing  various  combinations  of  powers 

2  2  5 

of  z  and  x.  We  are  interested  in  the  coefficient  of  the  terms  z  x  and  z  x. 

2  2  4 

Next  we  replace  z  with  z-u  and  pick  out  the  coefficients  of  u  x  and  uz  x.  We 


obtain  finally 
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«  4  2 

k202-Ta0D  (10a) 

kui  *  i5  a6°  D  »  tl0b) 

where  a  is  the  strength  parameter  of  the  AS  Morse  potential  and  D  is  the  well 

1  /2 

depth.  The  parameter  a  may  be  written  as  (uAS/3D)  •  To  the  best  of  our 

knowledge,  D  is  an  unknown  quantity  for  OH/SiC^,  but  we  shall  take  it  to  be  1 

-2 

eV,  noting  that  both  energy  and  phase  relaxation  depend  upon  D  .  Thus  changing 
D  by  a  factor  of  2  will  result  in  a  change  in  rate  of  a  factor  of  4 .  Finally, 
we  use  a  Debye  temperature  for  Si02  of  485  K  and  take  the  molecular  weight,  M, 

of  the  SiO.,  unit  cell  to  be  60.  From  Eqs.  (8a)  and  (9a)  we  now  obtain  T^  -  4ps 

4  -1  -4  -1 

and  T1  -  2  x  10  ps,  corresponding  to  widths  of  1.3  cm  and  2.6  x  10  cm  , 

respectively.  It  is  immediately  clear  that  the  calculated  value  for  T1  i3  about 

two  orders  of  magnitude  larger  than  the  experimental  value.1  However,  the  value 

of  T^  is  in  much  better  agreement  with  the  available  Information  on  this 

system,11  where  it  is  believed  that  the  pure  dephasing  width  <1  cm  1 . 

The  large  discrepancy  between  the  observed  and  calculated  values  of  T1  is 

very  interesting,  particularly  since  one  of  the  approximations  involved  the 

assumption  of  a  harmonic  adbond  for  the  calculation  of  k.^  in  Eq.  (9b). 

Removal  of  this  assumption  will  undoubtedly  increase  the  estimate  of  T1  still 

further.  It  is  unlikely  that  the  approximation  either  of  the  Debye  model  or  of 

the  linear  chain  dynamics  are  very  important  here.  Of  much  greater  importance 

is  the  crudity  of  our  use  of  a  Morse  potential  to  estimate  anharmonicities, 

particularly  as  (for  T^ )  the  relevant  terms  are  sextic. 

However,  the  discrepancy  is  sufficiently  large  that  we  may  tentatively 

postulate  a  different  mechanism  for  the  energy  relaxation.  The  reason  the 

mechanism  given  here  leads  to  such  long  relaxation  times  is  that  the  vibrational 
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quantum  of  the  molecular  bond  is  so  large  that  energy  transfer  requires  several 
quanta  of  excitation  of  other  modes,  and  this  is  a  process  requiring  at  least  a 
sextic  anharmonicity .  The  direct  excitation  of  the  substrate  requires  at  least 
eleven  phonons  and  may  be  presumed  to  be  insignificant.  It  is  possible  that  the 
OH  couples  to  another  kind  of  high-frequency  mode,  such  as  the  Si-OH  bend,  although 
this  mode  is  likely  to  be  of  a  somewhat  lower  frequency  than  the  Si-OH  stretch  we 
have  considered.  An  interesting  possibility  is  the  direct  energy  exchange  between 
adjacent  hydoxyl  groups  on  the  surface.  Even  for  weak  coupling,  such  a  pathway  may 
be  favored  because  the  vibrational  energy  mismatch  is  zero.  However,  the  irrevers¬ 
ibility  of  energy  relaxation  demands  several  adjacent,  anharmonically-coupled 
hydoxyl  groups.  In  this  regard,  therefore,  it  will  be  interesting  to  see  if 
increases  with  decreasing  coverage. 

As  far  as  pure  dephasing  is  concerned,  the  influence  of  an  adjacent  hydroxyl 

1  3 

may  be  crudely  modelled  in  terms  of  a  dipole-dipole  interaction.  Since  the 
dipole  moment  depends  upon  vibrational  state,  such  an  interaction  gives  rise  to 
a  relative  displacement  of  vibrational  levels,  which  leads  to  a  dephasing.  For 
H^O  the  vibrational  dependence  of  the  dipole  is  known  and  the  dephasing 
contribution  at  room  temperature  may  be  estimated  to  be  about  .01  cm  . 

Assuming  things  are  similar  for  OH,  we  can  therefore  rule  out  this  contribution 
to  the  pure  dephasing,  since  the  calculated  phonon-mediated  contribution  is  two 
orders  of  magnitude  faster.  The  calculation  of  the  contribution  of  this 
admolecule-admolecule  interaction  to  energy  relaxation  is  less  straightforward, 
and  will  be  possible  only  when  the  anharmonic  interaction  between  adjacent 
hydroxyl  groups  is  known.  It  is  also  possible  that  the  k^  term  is  much  larger 
than  our  estimate,  so  that  energy  exchange  with  the  lattice  is  important. 
Definitive  answers  await  further  experimental  information. 

During  preparation  of  this  manuscript  we  become  aware  of  a  recent  paper  by 
1 2 

Persson,  which  is  similar  in  spirit  to  the  present  work.  There  a  comparison 


was  made  between  energy  relaxation  and  dephasing  of  the  adbond,  specifically  for 
CO  on  nickel.  The  present  work,  on  the  other  hand,  describes  relaxation  of  the 
molecular  bond.  Moreover,  we  have  for  the  first  time  presented  explicit 
analytical  forms  for  T^. 

Conclusion 

We  have  derived  expressions  for  pure  vibrational  dephasing  and  for 
vibrational  energy  relaxation  of  a  diatomic  molecule  on  a  surface.  The  model 
involves  the  coupling  of  the  molecular  mode  to  bulk  phonons.  This  is  given  in 
terms  of  the  anharmonicities  of  the  molecular  and  admolecular  vibrations.  By 
making  an  estimate  of  thee  quantities  for  OH  on  SiO^,  we  have  shown  that  the 
prediction  of  the  pure  vibrational  phase  relaxation  is  in  agreement  with 
available  experimental  information,  but  the  energy  relaxation  rate  is  two  orders 
of  magnitude  smaller  than  experiment.  This  indicates  the  possibility  of 
another  mechanism  in  this  case,  and  the  obvious  candidate  is  energy  transfer 
between  adjacent  hydroxyl  groups.  A  definitive  answer  will  not  be  possible  until 
the  anharmonicities  are  known  accurately.  The  influence  of  surface  phonons  will 
be  assessed  at  a  later  t ime . 1 °  The  inclusion  of  such  modes  is  expected  to  be 
important  for  molecules  lying  across  the  surface  as  opposed  to  end-on,  as  in  the 
present  case.  This  is  because  of  the  expected  efficient  coupling  between  the 
lateral  motion  of  the  surface  atoms  and  the  motion  of  the  molecular  bond. 
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